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1.1 Production of antimatter in the galaxy 

The indirect detection of particle dark matter (DM) is based on the search 
for anomalous components in cosmic rays (CRs) due to the annihilation 
of DM pairs in the galactic halo, on the top of the standard astrophysical 
production. These additional exotic components are potentially detectable 
at Earth as spectral distortions for the various cosmic radiations: 

X + X^ qq^W^W" ^ p,D,e^ -fkv's . (1.1) 

Detection of the DM annihilation products has motivated the spectacular 
development of several new experimental techniques. They range from de- 
tectors on ballons or in space for the study of antimatter and gamma-rays, 
to large area cosmic-ray and gamma-ray detecors on the ground to neutrino 
telescopes underground for the study of the neutrino component. In the fol- 
lowing, we will discuss in detail the antimatter component of DM indirect 
searches, namely antiprotons, antideuterons, and positrons. 



1.2 Propagation of antinuclei in the Galaxy 

Whatever the mechanism responsible for their production, charged cosmic 
rays subsequently propagate through the galactic magnetic field and are 
deflected by its irregularities : the Alfven waves. In the regime where the 
magnetic turbulence is strong - which is the case for the Milky Way - Monte 
Carlo simulations |4qI) indicate that it is similar to space diffusion with a 
coefficient: 

K{E) = Kof3 {TZ/l GYf , (1.2) 

which increases as a power law with the rigidity TZ = pc/Ze of the parti- 
cle. In addition, because the scattering-centers drift inside the Milky Way 
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Fig. 1.1. Schematic edge-on view of the Milky Way diffusive halo (DH) as seen 
by a cosmic ray physicist. The stellar and gaseous disc is sandwiched between two 
thick layers which contain turbulent magnetic fields. After having been accelerated 
by SN driven shock waves or produced by DM species annihilating in the galactic 
halo, cosmic rays diffuse on magnetic inhomogeneities and are wiped away by a 
galactic wind with velocity Vc- They can lose energy and are also mildly subject to 
diffusive reacceleration. The former process is by far the dominant one in the case 
of electrons and positrons. This diagram has been borrowed from the review (jU). 



with a velocity 14 ~ 20 to 100 km , a second-order Fermi-mechanism 
is responsible for some mild diffusive reacceleration. Its coefficient Kee de- 
pends on the particle velocity (3 and total energy E and is related to the 
space diffusion coefficient K{E) through: 



Kee 



9 " 



K{E) 



(1.3) 
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Finally, galactic convection wipes cosmic rays away from the disc with a 
velocity Vc ~ 5 to 15 km s~^. We can assume steady state for the various 
populations of particles and write the mas ter e quation for the space and 
energy distribution function ip = dn/dE as (jl02l ): 

d,{Vcij) - i^A^ + dE{b'°'\E)^ - KEE{E)dE^] = q{^,E) . (1.4) 

This equation applies to any charged species - nuclei, protons, antiprotons 
or positrons - as long as the rates for production q and energy loss b^°^^{E) 
are properly accounted for. 

The solution of the master equation 11.41 has been deeply investigated and 
several different techniques lead to very similar fluxes at the Earth (f?^ : 



4J). One possibility is a completely numerical solution. This is the way 
followed in the Galprop model (j2ol : Is^ ) which, thanks to a realistic gas 
distribution, can calculated gamma rays in addition to charged species (jOSl ). 
A completely different approach of the CRs transport through the DH relies 
on the calculation of the Green function of Eq. 11.41 which describes the 
probability for a CR that is produced at a location x with energy Es to 
be detected at the Earth with a degraded energy E. The Green function is 
then integrated over the diffusive volume (DH) and the energy range. This 
is i.e. the approach chosen to propagate positrons (19; sl; 87). 

The Bessel expansion method, based on the cylindrical symmetry of the 
DH and on approximate values for the ISM (not relevant for charged CR 
propagation) permits a 2D fully analytical model. Numerical solution is 
re quir ed only for the diffusion in energy space. This is the model developed 
m I45I : 4ll ) and detailed more extensively in the following of this chap- 

ter. Very similar results for the propagation of stable primary and secondary 
nuclei have been obtained with the modified weighted slab technique (j27l ). 

According to this approach the region of the Galaxy inside which cosmic 
rays diffuse - the so-called diffusive halo (DH) - is pictured as a thick disc 
which matches the circular structure of the Milk Way as shown in Fig. 11.11 
The galactic disc of stars and gas, where primary cosmic rays are accelerated, 
lies in the middle. It extends radially 20 kpc from the center and has a half- 
thickness h of 100 pc. Confinement layers where cosmic rays are trapped 
by diffusion lie above and beneath this thin disc of gas. The intergalactic 
medium starts at the vertical boundaries z = ±L as well as beyond a radius 
of r = i? = 20 kpc. The half-thickness L of the diffusive halo is not known 
and reasonable values range from 1 to 15 kpc. The diffusion coefficient K is 
the same everywhere whereas the convective velocity is exclusively vertical 
with component Vc{z) = Vc sign(z). This galactic wind, which is produced 
by the bulk of disc stars like the Sun, drifts away from its progenitors along 
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the vertical directions. The normahzation coefficient Kq, the index (5, the 
galactic drift velocity Vc and the Alfven velocity Va may be determined by 
studying the boron-to-carbon ratio (B/C) which is quite sensitive to cosmic 
ray transport and which may be used efficiently as a constraint. 

The Bessel expansion method takes advantage of the axial symmetry of 
the DH and enforces a vanishing cosmic ray flux at a distance i? = 20 kpc 
from the rotation axis of the Galaxy. This condition is actually implemented 
naturally by the following series expansion for ^l^: 



The Bessel function of zeroth order Jq vanishes at the points ai . The radial 
dependence of ip is now taken into account by the set of its Bessel trans- 
forms Pi{z,E). The source term q may also be Bessel expanded into the 
corresponding functions Qi{z,E) so that the master equation (jl.4p becomes 



Here, energy loss and diffusive reacceleration are confined inside the galactic 
disc - which is considered infinitely thin, hence the presence of an effective 
term 2 h5{z). The form of the source terms Qi{z, E) which appear in equa- 
tion ()1.6p depends on the nature of the cosmic ray particle. 

1.3 Antiprotons in Cosmic Rays 

In the case of antiprotons, the following mechanisms can in principle con- 
tribute to the source term in the transport equation II. 4t 

• The spallation of high-energy primary nuclei impinging on the atoms 
of the interstellar medium inside the galactic disc produces secondary 
antiprotons. 

• The annihilation of DM candidate particles throughout the Milky Way 
halo generates primary antiprotons. Notice that WIMP annihilations take 
place all over the diffusive halo. 

• Tertiary antiprotons result from the inelastic and non-annihilating inter- 
actions with a nucleon at rest. The energy transfer may be sufficient 
to excite it as a A resonance. This mechanism redistributes antiprotons 
toward lower energies and flattens their spectrum (0). This yields the 




(1.5) 



1=1 
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source term: 

- apu^px{Ep) riH /3p ipp{r, Ep) , (1.7) 

where the inelastic and non-annihilating differential cross section in this 
expression can be approximated by: 

- ■ ^'-'^ 

The initial antiproton kinetic energy is denoted by Tp. In order to take 
into account elastic scatterings on helium, one simply has to replace the 
hydrogen density by nn + 4^/^ nne- 
• Antiprotons may also annihilate on interstellar H and He. This leads to a 
negative source term — ijj, where the annihilation rate Tp^^ is defined 
as 

■pann _ann n , ^ann o fl n\ 

The annihilation cross section <7p^ can be borrowed from (j^; 0) and mul- 
tiplied by a factor of 4-^/^ ~ 2.5, taking into account the higher geometric 
cross section, to get cr^^^. The average hydrogen nn and helium nue den- 
sities in the galactic disc are respectively set equal to 0.9 and 0.1 cm~^. 



1.3.1 Secondary antiprotons 

Secondary antiprotons are produced by CR rays-spallation on the inter- 
stellar medium. They represent the background when searching for small 
contributions coming from exotic sources, like signals from DM annihilation. 

The rate for the production of secondary antiprotons takes the following 
form: 

r+oo 

qf''{r,Ep) = / riH.a X /3p,„ Vp,a(?',^p,a) X f^£^p,a X 37r(^P,a -E^p) , 
Jeo „ dEp 



P,q: 



(1.10) 

for the interactions between cosmic ray protons and a particles, and hydro- 
gen and helium nuclei in the interstellar medium (ISM). 

Equation (jL^J may be solved according to the method outlined in Ap- 
pendix B of (|3ll ) for an antiproton source located only in the galactic disk, 
as is the case for secondary antiprotons. 

Some approximations may be at hand: in particular, setting the energy 
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Fig. 1.2. Left panel: IS antiproton flux for the B/C best flt model and two param- 
eterizations of the production cross section. Right panel: propagation uncertainty 
envelopes of the IS p/p ratio for the same production cross sections as in the left 
panel. All data are demodulated using the force-field approximation: AMS 98 (37). 
IMAX 92 Jll), CAPRICE 94 (0), WIZARD-MASS 91 dH), CAPRICE 98 ®, 
BESS93 (H^BESS 95+97 (28^, BESS 98 (33), BESS 99 and 2000 (38), BESS 2002 
dH), BESS Polar ([zi), WIZARD-MASS 1 Q^), HEAT-p ((2^), and PAMELA @. 



loss rate 6^°^^ and the energy diffusion coefficient Kee equal to zero does 
not affect sizeably tlie solution of the diffusion equation. 

As for the the cosmic ray proton and helium fluxes, they can be borrowed 



from (j3ll ). where a fit to high energy (> 20 GeV/n) data is proposed. A 



more recent parameterization on IS fluxes has been derived in (19 ll ). 

Once obtained the interstellar (IS) fluxes of antiprotons at the Sun's posi- 
tion in the Galaxy, we have to further propagate them inside the heliosphere, 
where the cosmic ray particles which eventually reach the Earth are affected 
by the presence of the solar wind. We model the effect of solar modulation 
by adopting the force field approximation of the full transport equation (0). 
In this model, the top-of-atmosphere (TOA) flux for a cosmic species <I>"'"'-'^ 
is obtained as: 



$IS(^IS) I pis 



(1.11) 



where E and p denote the total energies and momenta of interstellar and 
TOA antiprotons, which are related by the energy shift: 

^TOA = ^IS _ ^ (1 12) 

where the parameter (p is determined by fits on cosmic ray data. A value 
(j) = 500 MV is indicative of periods of minimal solar activity. 

The secondary IS p flux is displayed in the left panel of Fig. 11.21 along 
with the data demodulated according to the force-field prescription. We 
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Halo model 


a 




7 


Ps [10^ Mq kpc-^] 


rs [kpc] 


Cored isothermal (5) 


2 


2 





7.90 


4 


NFW 97 (18) 


1 


3 


1 


5.38 


21.75 


Moore 04 (50) 


1 


3 


1.16 


2.54 


32.62 



Table 1.1. Parameters in equation il.l5\) for different halo models. The 
scale radius rs and density ps are strongly correlated with the virial mass of 
the Galaxy l32) and the values are borrowed from iSdi) for the Milky Way. 

When the DM distribution is cuspy (7 > 1 j the divergence at the galactic 
center is smoothed according to the prescription of ifyj). 

either use the DTUNUC (l3ll ) p production cross sections (solid line) or 



those discussed in Refs. (|6ll : l75l ) (dashed line). The differences between 
the two curves illustrate the uncertainty related to the production cross 
sections, as emphasized in (jsil). where a careful and conservative analysis 
within the DTUNUC simulation settled a nuclear uncertainty of ~ 25% over 
the energy range 0.1 — 100 GeV. The conclusion is similar here, although 
the two sets of cross sections differ mostly at low energy. In the right panel, 
along with the demodulated p/p data, we show the curves bounding the 
propagation uncertainty on the p calculation based either on the DTUNUC 



(j3ll ) p production cross sections (solid lines) or those borrowed from (j75l ) 



(dashed lines). The uncertainty arising from propagation is comparable to 



the nuclear one (|3ll). The detailed calculation of that secondary component 
(0) has required the determination of thepropagation-diffusion parameters 
that are consistent with the B/C data (jSSi ). By varying those parameters 
over the entire range allowed by the cosmic~ray nuclei measurements, the 
theoretical uncertainty on the antiproton secondary flux has been found to 
be 9% from 100 MeV to 1 GeV. It reaches a maximum of 24% at 10 GeV 
and decreases to 10% at 100 GeV. 

From Fig. 11.21 it is manifest that the secondary contribution alone ex- 
plains experimental data on the whole energetic range. It is not necessary 
to invoke an additional component to the standard astrophysical one. 



1.3.2 Antiprotons from DM annihilation 

The antiproton signal from annihilating DM particles leads to a primary 
component directly produced throughout the DH. The differential rate of 
production per unit volume and time is a function of space coordinates and 
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Case 


5 


A'o [kpcVMyr] 


L [kpc] 


Vc [km/s] 


Va [km/s] 


MIN 


0.85 


0.0016 


1 


13.5 


22.4 


MED 


0.70 


0.0112 


4 


12 


52.9 


MAX 


0.46 


0.0765 


15 


5 


117.6 



Table 1.2. Typical combinations of diffusion parameters that are compatible 
with the B/C analysis fs^)- As shown in i5i). these propagation models 
correspond respectively to minimal, medium and maximal primary 

antiproton fluxes. 



antiproton kinetic energy Tp. It is defined as: 

gD^(r,z,r,) = (^^) ' (1-13) 

where (crann'w) denotes the average over the Galactic velocity distribution 
function of the WIMP pair annihilation cross section o"ann multiplied by the 
relative velocity v. For a relic CDM particle able to explain the observed 



amount of cosmological dark matter (jOSl : 1971 : l93l ) its value falls in the range 
2-^3- 10~^^ cm^ s~^ (unless special situations occurs, like dominant p-wave 
annihilation, dominant coannihilations or modified cosmology). For any DM 
candidate, o"ann is nevertheless calculated from the model parameters. The 
WIMP mass is denoted by cJann and p^if, z) is the mass distribution function 
of DM particles inside the Galactic halo. The quantity ^ parameterizes the 
fact that the dark halo may not be totally made of the species under scrutiny 
when this candidate possesses a relic abundance which does not allow it to be 
the dominant DM component (see e.g (jssl ) or (fiol )). In this case ^ < 1. The 
quantity g{Tp) in Eq. (jl.l3p denotes the antiproton differential spectrum per 
annihilation event, defined as: 

^ = E<^^- (1-14) 

The annihilation into a quark or a gluon h is realized through the various 
final states F with branching ratios -B^^ • Quarks or gluons may in fact be 
directly produced when a WIMP pair annihilates or they may alternatively 
result from the intermediate production of Higgs bosons or gauge bosons. 
Each quark or gluon h then generates jets whose subsequent fragmentation 
and hadronization yield an antiproton energy spectrum dNp/dEp. The sin- 
gle production spectra are usually evaluated within Monte Carlo simulations 
of electroweak annihilation events (js^ ) . For details on the calculation of pure 
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TI" (GeV) 

Fig. 1.3. Top-of-atmosphere antiproton fluxes as a function of the antiproton ki- 
netic energy for the = 100 GeV reference case. The upper (lower) set of curves 
refer to the maximal (minimal) set of of astrophysical parameters. Solid curves 
show the interstellar fluxes. Broken curves show the effect of solar modulation at 
different periods of solar activity: (p = 500 MV (long dashed), (p — 700 MV (short 
dashed), = 1300 MV (dotted). 



state dNp/dEp and of g{Tp) in a Minimal Supersymmetric extension of the 
Standard Model (MSSM), see the Appendix of Ref. <^). 

The distribution of DM inside galaxies is a very debated issue. Different 
analyses of rotational curves observed for several types of galaxies strongly 
favour a cored dark matter distribution, flattened towards the central regions 
(Ref. 0) and references therein). On the other side, many collisionless 
cosmological N-bod y si mulations in A-CDM models are in good agreement 
among themselves (|9J), but for the very central regions some resolution 
issues remain open. It has been recently stressed that asymptotic slopes may 
not be reached at all at small scales (57; 7^; 71; 66; 76). However, it is not 
clear whether the central cusp is steepened or flattened when the baryonic 
distribution is taken int o ac count (e.g . (jgl; 70)). A detailed discussion 
may be found in Ref. ([lOl|; S Eool ). We wish to stress that the DM 



distribution is a crucial ingredient for the indirect signal into gamma-rays 
or neutrinos, while for antimatter it has been shown to be less relevant 
(|5ll : |89j) , since antimatter diffuses and therefore the signal at Earth is more 
local, less dependent on the inner structure of the Galaxy as compared to 
gamma-rays and neutrinos. 

To be definite, we will consider here a spherical DM galactic distribution 
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Fig. 1.4. Left: The scatter plot shows the antiproton flux at solar minimum from 
neutralino annihilation calculated at Tp = 0.23 GeV, as a function of the neu- 
tralino mass for a generic scan of in a low-energy MSSM and for the MED set of 
astrophysical parameters (j52il . Crosses (in red) refer to cosmologically dominant 
neutralinos (0.05 < fi^/i^ < 0.3); dots (in blue) refer to subdominant relic neutrali- 
nos (ri/i^ < 0.05). Right: The same as in the left panel, but calculated for in scan 
of supersymmetric framework where gaugino non-universality is not assumed, and 
therefore lighter neutralinos are present . The astrophysical parameters are set 
at the MIN case, and solar modulation is at solar minimum. 



with a dependence on galactocentric distance r parameterized by: 

(7-/3)/" 



p{r) 



Ps 



1 + 



(1.15) 



The different profiles of Table 11.11 basically span the whole range of reason- 
able halo models with respect to indirect dark matter detection prospects. 

The solution of the transport equation for a generic source term qP"'^(r, z, E) 
within Bessel expansion has been derived in (j39j): 



A^f'P""^(z) 



exp 



Vc{\z\-L) 
2K 



y^{L) 



Ai sinh(5iL/2) 



cosh(5,z/2) + ' / sinh(5,z/2) 



where: 



exp 



(^^" - "'0 """"^ (t " ^ (1.17) 
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and the quantities Si and Ai are defined as: 

f T/2 ^2 ] 1/2 ( Q T ^ 

5i = |-^ + 4|2 and Ai{E)^2hTf^ + V, + KS,coth^-^j 

In particular, at z = where fluxes are measured, we have: 

A'r""(0)=exp(^)^-M|^^ (1.19) 
The three propagation models featured in Table [L2] have been drawn from 



(1.18) 



(jSlI ). The MED configuration provides the best fit to the B/C measurements 
whereas the MIN and MAX models point respectively to the minimal and 
maximal allowed antiproton fluxes which can be produced by WIMP anni- 
hilation. 

Fig. 11.31 shows the TOA antiproton fluxes for a reference WIMP with 
= 100 GeV, 2.3 • 10~^^ cm^ s"^, pure bb annihilation final state and for 
the maximal and minimal sets of astrophysical parameters. The variation 
of the astrophysical parameters induces a much larger uncertainty on the 
primary than on the secondary flux: in the first case, the uncertainty reaches 
two orders of magnitude for energies Tp < 1 GeV, while in the second case 
it never exceeds 25%. The figure shows that solar modulation has the effect 
of depleting the low-energy tail of the antiproton fiux. The effect is clearly 
more pronounced for periods of high solar activity, when the solar wind is 
stronger. 

Predictions for antiprotons in various realization of sup er symmetric the- 
ories have been performed (see. e.g. (52; 54; 4^; 60; 65; 0) 

. In Fig. I1.4l we show the predicted antiproton fiux at kinetic energiy Tp = 
0.23 GeV, for two different supersymmetric models: the left panel (jM) refers 
to a low-energy realization of the MSSM, while the right panel (j52l ) stands 
for a scan of a supersymmetric scheme where gaugino non-universality is not 
assumed and therefore light neutralinos are present (|48l : 1471) . In both cases 
the shaded (yellow) region denotes the amount of antiprotons, in excess 
of the secondary component , which can be accommodated at Tp = 0.23 
GeV in order not to exceed the observed flux, as measured by BESS. All 
the points of the scatter plot that lie below the horizontal black line are 
therefore compatible with observations. 



1.4 Antideuterons in Cosmic Rays 

In the seminal Paper it was proposed to look for cosmic antideuterons 
(D) as a possible indirect signature for galactic dark matter. It was shown 
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that the antideuteron spectra deriving from DM annihilation is expected 
to be much flatter than the standard astrophysical component at low ki- 
netic energies, < 2-3 GeV/n. This argument motivated the proposal of 



a new space-borne experiment (j46|; |56|; l67l ) looking for cosmic antimatter 
(antiproton and antideuteron) and having the potential to discriminate be- 
tween standard and exotic components for a wide range of DM models. The 
present experimental upper limit ()62l ) is still far from the expectations on 
the secondary antideuteron flux which are produced by spallation of cosmic 
rays on the interstellar medium (1^; 61)) but perspectives for the near future 
are very encouraging. 

Once the astrophysical framework for the transport of (anti)nuclei is set, 
the calculation of the antideuteron flux rests on the D specificities regarding 
the source term (i.e. primary or secondary) and its nuclear interactions 
[r.h.s. of Eq. ()1.4p ]. The tertiary term (see Eq. II. 7p requires specific inelastic 
non-annihilating cross sections which are detailed in the Appendix of Ref. 

(B). 

The production of cosmic antideuterons is based on the fusion process of 
a p and h pair. One of the simplest but powerful treatment of the fusion 
of two or more nucleons is based on the so-called coalescence model which, 
despite its simplicity, is able to reproduce remarkably well the available data 
on light nuclei and antinuclei production in different kinds of collisions. In 
the coalescence model, the momentum distribution of the (anti)deuteron is 
proportional to the product of the (anti)proton and (anti)neutron momen- 
tum distribution (jl|;y)- That function depends on the difference between 
(anti)nucleon momenta. It is strongly peaked around ~ (compare the 
minimum energy to form a D, i.e. 4mp, with the binding energy ~ 2.2 MeV), 
so that 

h 

2 



The D density in momentum space is thus written as the p density times 
the probability to find an n within a sphere of radius po around kp (see, e.g. 
Ref. (nl)): 

dkj~, 3 d}^- (ikfi 

The coalescence momentum po is a free parameter constrained by data on 
hadronic production In (61), a large set of data is used, including many pA 
reactions (see their Tab. I and references therein). This leads to an estimate 
of po = 79 MeV. At LEP energies, (anti)deuteron production occurs through 
e'^e" annihilations into qq pairs, an electroweak mechanism similar to the D 
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Fig. 1.5. Left panel: Contribution of all nuclear channels to the D secondary flux. 
Dashed lines, from top to bottom refer to: p+H, p+He, He+H, He+He. Dotted 
lines, from top to bottom stand for: p+H, p+He. Solid line: sum of all the compo- 
nents. Right panel: Dominant uncertainties on the interstellar secondary D flux. 
Solid lines: propagation uncertainty band. Dotted lines: nuclear uncertainty band. 



production in DM annihilation reactions. Based on theoretical arguments, it 
has been argued (jlll ) that the antideuteron yields in e^e~ reactions should 
be smaller than in hadronic reactions. However, the ALEPH Collaboration 
(j72l ) has found that this theoretical prediction (see Fig. 5 in ALEPH paper) 
underestimates their measured D inclusive cross section. They derive (see 
their Fig. 6) a value B2 = 3.3 ± 0.4 ± 0.1 x 10"^ GeV^ at the Z resonance, 
which translates into po = 71.8 ± 3.6 MeV, very close to the po = 79 MeV 
derived for the hadronic production. From a rough estimate of the fusion 
process, one can estimate that the antideuteron fluxes are a factor of 10^ 
lower than the antiproton ones. We will see in the following, and by com- 
paring the results of the previous Section, that this result is recovered for 
both the primary and the secondary D fluxes at Earth. 



1.4-1 Secondary Antideuterons 

The secondary D flux is the sum of the six contributions corresponding to p, 
He and p cosmic ray fluxes impinging on H and He IS gas (other reactions 
are negligible ()6ll )). Contributions to the D flux from p + H and p + He 
reactions are evaluated using the p flux calculated in the same run. The 
pro duction cross sections for these specific processes are those given in Ref. 
(j6ll ). The solution to the propagation equation has the same expression as 
for secondary antiprotons (j92l ). 

The different contributions to the total secondary antideuteron flux, cal- 
culated for the best fit propagation configuration (the "MED" one in Table 
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ll.2p . are shown in the left panel of Fig. 11.51 As expected, the dominant 
production channel is the one from p-H collisions, followed by the one from 



cosmic protons on IS helium (p-He). As shown in Ref. (j6ll ). the p+H chan- 
nel is dominant at low energies, and negligible beyond a few GeV/n. The 
effect of energy losses, reacceleration and tertiaries add up to replenish the 
low energy tail. The maximum of the total flux reaches the value of 2 • 10^^ 
particles (m^ s sr GeV/n)~^ at 3-4 GeV/n. At 100 MeV/n it is decreased by 
an order of magnitude, thus preserving an interesting window for possible 
exotic contributions characterized by a flatter spectrum (see next Section). 

For the determination of the propagation uncertainties, one can calculate 
the secondary antideuteron flux for all the propagation parameter combina- 
tions providing an acceptable fit to stable nuclei (|35l). The resulting envelope 
for the secondary antideuteron flux is presented in Fig. 1 1.51 - left. The solid 
lines delimit the uncertainty band due to the degeneracy of the propaga- 
tion parameters: at energies below 1~2 GeV/n, the uncertainty is 40-50 % 
around the average flux, while at 10 GeV/n it decreases to ~ 15 %. This 
behavior is analogous to that obtained for p (jsil ) and is easily understood. 
The degenerate transport parameters combine to give the same grammage 
in order to reproduce the B/C ratio. Indeed, the grammage crossed by C 
to produce the secondary species B is also crossed by p and He to produce 
the secondary p and D. In short, a similar propagation history associated 
with a well constrained B/C ratio explains the small uncertainty. The pos- 
sible nuclear uncertainty can arise from two different sources. The first one 
is directly related to the elementary production process da^ and has been 
estimated to be conservatively ±50% (jool ). Second, there is the uncertainty 
on the coalescence momentum pQ. Using an independent model (i.e. differ- 



ent from the coalescence scheme) for D production, Ref. (j6ll ) found that, 
conservatively, the D background was certainly no more than twice the flux 
calculated with po = 79 MeV. The dotted lines in Fig. II. 5t right panel, 
take into account the sum of all the possible uncertainties of nuclear source 
(jool ). At the lowest energies the flux is uncertain by almost one order of 
magnitude, at 100 GeV/n by a factor of 4. 



1.4-2 Antideuterons from DM annihilation 

The source term for primary D has the same form as the one for primary 
antiprotons in Eq. 11.131 [p ^ D). The production of antideuterons from the 
pair-annihilation of supersymmetric dark matter particles in the halo of our 
Galaxy was proposed in |26l ). and subsequently discussed in (jsol) also for 
universal extra-dimension, Kaluza-Klein and warped extra-dimension dark 
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matter models. In ref. (]90|) primary antideueterons in different supersym- 
metric scenarios have been calculated in the full 2D propagation model, with 
a thorough estimation of all the possible uncertanty sources. As previously 
discussed (see Sect. II. 4p the production of a D relies on the availability of 
a p ~ h pair in a single DM annihilation. As in Ref. (0), we assume that 
the probability to form an antiproton (or an antineutron) with momentum 
kp {kn)i is essentially isotropic: 

^(X + X^P+-..) = ^TikpEpFp{V~s = '2m,Ep) . (1.22) 

Applying the factorization-coalescence scheme discussed above leads to the 
antideuteron differential multiplicity 

F,h 




We assume, as discussed in Sect. 11.41 that the same value of the coalescence 
momentum pQ = 79 MeV holds as for hadronic reactions. 

The solution to the diffusion equation is provided by Eq. I1.19j where the 
specific cross sections have been reviewd, for instance, in Ref~(9(?). In the 
left panel of Fig. 11.61 the secondary D flux for the median configuration of 
Table [L2] is plotted alongside the primary flux from =50 GeV, calculated 
for the maximal, median and minimal propagation scenarios. The present 
BESS upper limit on the (negative) antideuteron search Q is at a level of 
2-10"'^ (m^ s sr GeV/n)"^. We also plot the estimated sensitivities of the 
gaseous antiparticle spectrometer GAPS on a long duration balloon flight 



(LDB) and an ultra-long duration balloon mission (ULDB) (|56l : l67l : l96l ). 
and of AMS-02 for three years of data taking. The perspectives to explore a 
part of the region where DM annihilation are mostly expected (i.e. the low- 
energy tail) are very promising. If one of these experiments will measure at 
least 1 antideuteron, it will be a clear signal of an exotic contribution to the 
cosmic antideuterons. 

The discrimination power between primary and secondary D flux may be 
deduced from the right panel of Fig. 11.61 The ratio of the primary to total 
TOA D flux is plotted as a function of the kinetic energy per nucleon, for 
the three representative propagation models and different WIMP masses 
(the annihilation cross section is again fixed at the reference value). This 
ratio keeps higher than 0.7 for < 1 GeV/n except for =500 GeV. For 
propagation models with L ^ 4 kpc - which is a very reasonable expectation 
- this ratio is at least 0.9 for masses below 100 GeV. Increasing the WIMP 
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Tj (GeV/n) T5 (CeV/n) 

Fig. 1.6. Left panel: TOA primary (red solid lines) and secondary (black dashed 
line) antideuteron fluxes, modulated at solar minimum. The signal is derived for 
a m^=50 GeV WIMP and for the three propagation models of Table 11.21 The 
secondary flux is shown for the median propagation model. The upper dashed 
horizontal line shows the current BESS upper limit on the search for cosmic an- 
tideuterons. The three horizontal solid (blue) lines are the estimated sensitivities for 
(from top to bottom): AMS-02^ GAPS on a long (LDB) and uhra-long (ULDB) 
duration balloon flights ((56l : 1671 : |9^. Right panel: Ratio of the primary to to- 
tal (signal-|-background) TOA antideuteron flux. Solid (black) curve refers to a 
WIMP mass of m^=50 GeV and for the MED propagation parameters. Dotted 
(black) lines show the MAX (upper) and MIN (lower) cases. Dashed lines refer 
to the MED propagation parameters and different masses, which are (from top to 
bottom): m^—lO, 100, 500 GeV (red, blue, magenta respectively). 

mass, we must descend to lower energies in order to maximize the primary- 
to-secondary ratio. However, for a =500 GeV WIMP we still have a 
50-60% of DM contribution in the 0.1-0.5 GeV/n range. Of course, the 
evaluation of the theoretical uncertainties presented in this Paper must be 
kept in mind while confronting to real data. Fig. 11.61 clearly states that the 
antideuteron indirect DM detection technique is probably the most powerful 
one for low and intermediate WIMP-mass haloes. 

The uncertainties due to propagation are similar to the ones quoted for 
primary antiprotons (jool ). At the lowest energies of hundreds of MeV/n the 
total uncertainty reaches almost 2 orders of magnitude, while at energies 
above 1 GeV/n it is about a factor of 30. 

Examples for theoretical calculations for neutralino dark matter are shown 
in Fig. II. 7| where the reaching capabilities of the GAPS detector are derived 
and shown in terms of the neutralino annihilation parameters. The shaded 
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Fig. 1.7. Left: GAPS Ultra Long Duration Flight (ULDF) reach compared to 
predictions for neutralino dark matter in low-energy supersymmetric models, shown 
in the plane effective annihilation cross section £,'^{<Ja,nnv)o vs. neutralino mass 
(jOOl). The solid, long-dashed and short-dashed lines show our estimate for the 
capability of GAPS ULDB of measuring 1, 10 and 100 events, respectively, for 
the MED propagation model. The scatter plot reports the quantity ^^(crann^^)o 
calculated in a low-energy MSSM (for masses above the vertical [green] dashed 
line) and in non-universal gaugino models which predict low-mass neutralinos (j47l : 
HI). [Red] Crosses refer to cosmologically dominant neutralinos, while [blue] dots 
stand for subdominant neutralinos. Grey point are excluded by antiproton searches. 
Right: The same as in the left panel, but for SUGRA scheme with non-universality 
in the higgs sector (pol ) . 



band denotes the region which will be covered by GAPS with a Ultra Long 
Duration Flight: sensitivity to neutralinos from a mass of a few GeV up 
to few hundreds of GeV is at hand and it will cover a large portion of the 
supersymmetric parameter space for models with gaugino non-universality 



(147|; |48|) or low-energy MSSM (jgol) 



1.5 Positrons in Cosmic Rays 

In the case of positrons and electrons, the master equation (jl.4p describing 
the propagation of cosmic rays throughout the DH is dominated by space 
diffusion and energy losses. Above a few GeV, synchrotron radiation in the 
galactic magnetic fields as well as inverse Compton scattering on stellar light 
and on CMB photons dominate, hence the positron loss rate 

^ioss(^) = /^\ = - . (1.24) 
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The energy of reference Eq is set equal to 1 GeV while the typical energy 
loss time te is of order 10^^ s. The master equation for positron propagation 
simplifies into 

- K^i, + dE[b^'''\E)i,] = q{^,E) . (1.25) 

Above a few MeV, positrons are ultra-relativistic and the rigidity TZ is pro- 
portional to the energy E. The space diffusion coefficient boils down to 
K{e) = Koe^ where e = E/Eq. The solution of equation (|1.25p proposed 



by (13; I21I ) is based on replacing the energy E by the pseudo-time 

i{E) = TE i^v{E) = ^1 , (1.26) 
and leads to the well-known heat equation 

^ - KoAi; = g(x,t) , (1.27) 

where the space and energy positron density is now given hy ip = if) 
whereas the positron production rate has become q = e^~^ q. Because e is 

dimensionless, both ■0 and q have the same dimensions as before. In this 
formalism, the energy losses which positrons experience are described by 
an evolution of their density with respect to the pseudo-time i. Equipped 
with these notations, we can write the cosmic ray positron density as the 
convolution 



i;,+{K,E)= / dEs d^'^s G,+{x,E ^xs,Es) qM^s^Es) . (1.28) 

JEa=E J DH 



rEs-- 

lEs=E J DH 

The positron propagator Gf,+{x,E <— xs,Es) is defined as the probability 
for a particle injected at X5 with the energy Es to reach the location x with 
the degraded energy E < Es- It is proportional to the Green function G of 
the heat equation (|1.27p through 

Ge+(x, E ^ xs, Es) = G (x, i ^ X5, is) , (1.29) 

where the connection between the energy E and pseudo-time t is given by 
relation ([TJU]). 

The build up our intuition, it is useful to derive the heat Green function 
G connecting a source at xg to the Earth in the simple case of an infinite 
DH. In this limit, there are no vertical boundaries and we get the Gaussian 
distribution 

6(xe,f^x.,f.) = {j-i^}"%xp|-j|i.| . (1.30) 
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where f = t — ts is the typical duration over which the positron energy 
decreases from Es to E. That timescale also includes information on the 
diffusion process. The distance between the source at and the Earth is 
denoted by r^. The concept of positron horizon is based on the Gaussian 
distribution ()1.30p which is roughly constant within a sphere of radius 



and decreases sharply outside. So does the positron Green function 0^+. 
The positron sphere - whose center is at the Earth where the observer stands 
- delineates actually the region of the diffusive halo from which positrons 
predominantly originate. The typical diffusion length Ad gauges how far 
particles produced at the energy Es travel before being detected with the 
energy E. It encodes at the same time the energy loss process and the 
diffusion throughout the magnetic fields of the Galaxy. A rapid inspection of 
equation (|1.26p shows that Ad increases as the detected energy E decreases, 
except for energies Es at the source very close to E. The positron sphere 
is indeed fairly small at high energies, say above ~ 100 GeV, whereas it 
spreads over several kiloparsecs below 10 GeV. In the case where Es = 100 
GeV for instance, Ad exceeds 3 kpc below an energy E' of ~ 8 GeV. 

The diffusive halo inside which cosmic rays propagate before escaping into 
the intergalactic medium is actually finite and G should account for that ef- 
fect. In spite of the boundaries at r = i? = 20 kpc, we can decide that 
cosmic ray diffusion is not limited along the radial direction and that it op- 
erates as if it took place inside an infinite horizontal slab with half-thickness 
L. Sources located beyond the galactic radius R are disregarded since the 
convolution ()1.28p is performed only over the DH. Because their energy is 
rapidly degraded as they propagate, positrons are produced close to where 
they are observed. Neglecting the effect of radial boundaries on the prop- 
agator Gg+ turns out to be a fair approximation (jssi ) because positrons do 
not originate from far away on average. The effects of the radial boundaries 
down at the Earth are not significant insofar as cosmic rays tend to leak 
above and beneath the diffusive halo at 2: = itL instead of traveling a long 
distance along the galactic plane. 

The infinite slab hypothesis allows the radial and vertical directions to be 
disentangled and to express the reduced propagator G as 



Ad = V4:Kof 



(1.31) 



G{x,i ^xs,is) 



1 



exp 



( 



V{z,t^zs,is) , (1.32) 



where the radial distance between the source at X5 and the point x of 
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observation is now defined as 

r = \{x - xsf + {y ~ ysf^ ^ ■ (1-33) 

Should the half-thickness L be very large, expression (jl.32p would boil down 
to the Gaussian distribution (jl.30p and the vertical propagator V would be 
given by the ID solution Vid of the heat equation (|1.27|) 



V {z,t^ zs,is) = ViD {z,i^ zs,is) = \, - exp \ - J^'^] 



(1.34) 



But the diffusive halo has a finite vertical extent. We need to implement 
the corresponding boundary conditions and impose that the positron density 
vanishes at z = ±L. 

(i) A first approach relies on the method of the so-called electrical images 
and has been discussed in (l2ll ). Any point-like source inside the slab is 
associated to the infinite series of its multiple images through the boundaries 
at z = ibL which act as mirrors. The n-th image is located at 

Zn = 2Ln + (-1)" zs , (1.35) 

and has a positive or negative contribution depending on whether n is an 
even or odd number. When the diffusion time f is small, the ID solu- 
tion (|1.34|) is a quite good approximation. The relevant parameter is actu- 
ally 

Jl_ = Jl 

AKof ~ Ad' 

and in the regime where it is much larger than 1, the propagation is insen- 
sitive to the vertical boundaries. On the contrary, when C, is much smaller 
than 1 , a large number of images need to be taken into account in the sum 

+ 00 

V {z,i ^ zs.is) = E (-1)" ViD(^,t ^2n,ts) , (1.37) 

n = — oo 

and convergence may be a problem. 

(ii) It is fortunate that a quite different approach is possible in that case. 
The ID diffusion equation ()1.27p actually looks like the Schrodinger equation 
- though in imaginary time - that accounts for the behaviour of a particle 
inside an infinitely deep ID potential well which extends from z = —L to 
z = +L. The eigenfunctions of the associated Hamiltonian are both even 

^n{z) = siYi{kn{L-\z\)} (1.38) 



C = = — , (1-36) 
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and odd 

ifUz) = sin {k'„iL-z)} (1.39) 

functions of the vertical coordinate z. The wave- vectors kn and k'^ are 
respectively defined as 

/ 1 \ TT TT 

kn = ^ ~ o ) T (even) and k'^ = n— (odd) . (1-40) 

\ Z J Li Li 

The vertical propagator may be expanded as the series 

V [zj *^ zsM) = ^ i |e~-^"^V'nfe) V'nC^) + e"^"^</:'U^s) , (1-41) 

n = 1 

where the time constants A„ and are respectively equal to i^o ^n' and 
Kq k'^ . In the regime where is much smaller than 1, for very large values 
of the diffusion time f, just a few eigenfunctions need to be considered in 
order for the sum ()1.4ip to converge. A close examination of these various 
expressions for G indicate that the energies E and Es always come into 
play through the diffusion length Ad so that the positron propagator may 
be written as 

G,+(x, E ^ X5, Es) = G (x ^ xs; Ad) . (1.42) 



1.5.1 Secondary Positrons 

Like for antiprotons, a background of secondary positrons is produced by the 
spallation of the interstellar medium by impinging high-energy particles. In 
that respect, the Milky Way looks like a giant accelerator where cosmic rays 
play the role of the beam whereas the galactic disc and its gas behave as the 
target. The dominant mechanism is the collision of protons with hydrogen 
atoms at rest producing charged pions vr^ which decay into muons /i^. 
The latter are also unstable and eventually lead to electrons and positrons 
through the chain 

p + H ^ X + 7r± 

fi^ > + Ve + . 

Below ~ 3 GeV, one of the protons is predominantly excited to a A res- 
onance which subsequently decays into either a neutral or a charged pion. 
The former species produces gamma rays whereas the latter particle de- 
cays into positrons. Above ~ 7 GeV, pion production is well described in 
the framework of the scaling model. Various parameterizations are given in 
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the hterature (j3;@) for the Lorentz invariant (LI) cross section E-j^ cfa/d^pj^. 
Positrons may also be produced through kaons although this channel is rare. 
The positron production cross sections of these processes have been carefully 
computed in the appendices of Notice also that useful parametric ex- 
pressions for the yield and spectra of the stable secondary species produced 
in p-p collisions have been derived from experimental data and summarized 
in (j68l ). 

T. Delahaye et al, (2008) 
= 1 I ' I ' I ' 

AMS 98 CAPRICE 94 HEAT 94 + 95 MASS 91 



IS 




Positron energy E^, [GeV] 



Fig. 1.8. Secondary positron flux as a function of the positron energy. The blue 
hatched band corresponds to the cosmic ray propagation uncertainty on the inter- 
stellar prediction whereas the yellow strip refers to the top-of-atmosphere fluxes. 
The long-dashed curves feature our reference model MED with a differential pro- 
duction cross section borrowed from (|68l ) and the most recent measurements by 
BESS (7?) of the cosmic ray proton and helium fluxes. The MIN, MED and 
MAX propagation parameters are displayed in table 11.21 The observations by 
CAPRICE ^ (green squares), HEAT ^ (blue triangles), MASS gl) (open 
circles) and AMS (|24l : [tI) (red dots) are also indicated. 

Cosmic ray protons with energy Ep induce a production of positrons per 
hydrogen atom with a rate 

dT'^{Ee) = ^{Ep ^ Ee) x/3pX {dnp = Vp(^p) x dEp} . (1.43) 
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This leads to the source term 

(x, ) = 4 vr nj^(x) J $p(x, Ep ) x dEp x ^{Ep ^ E,) . (1.44) 

That relation can be generalized in order to incorporate cosmic ray helium 
nuclei as well as interstellar helium. The gas of the galactic plane is generally 
assumed to be homogeneously spread. Because positrons detected at the 
Earth originate mostly from the solar neighborhood, we can safely disregard 
the space dependance of the proton and helium fluxes. Making use then of 



the measurements (13l|; 1771 ) in relation (|1.44p , we can express the secondary 
positron flux as 

<^l%%Q,e^EjEo) = ^x^x r^desxi{Xj,)xqlf{Q,es) . (1.45) 

4 7r Je 

The integral / is the convolution of the reduced positron Green function G 
over the galactic disc alone 

/(Ad) = / (i'x5 G(x0^X5; Ad) , (1.46) 

J disc 

and depends only on Ad as discussed above. 

This method has been recently used js^ ) to derive the positron flux fea- 
tured in figure 11.81 At 1 GeV, the width of the IS uncertainty strip corre- 
sponds to an increase by a factor of ~ 6 between the smallest and the largest 
positron fluxes allowed by the B/C constraint. That factor decreases down 
to 3.9 at 10 GeV and reaches 2.9 at 100 GeV. Once modulated with a Fisk 
potential cpp oi 600 MV, the blue hatched region is transformed into the 
yellow TOA band. Quite surprisingly, the MIN model (red solid curve) cor- 
responds now to the maximal secondary positron flux whereas the MAX 
configuration (red short-dashed line) yields the minimal prediction. At 
fixed energies e and es, the smaller the diffusion coefficient Kq, the smaller 
the diffusion length Ad and the larger the integral ()1.46p . The latter reaches 
its maximal value of 1 whenever Ad is smaller than the disc half-thickness 
h of 100 pc. 



1.5.2 DM signals in Cosmic Positrons 

An excess of the positron measurements with respect to the astrophysical 
background seems to appear above 10 GeV. This trend was present in the 
HEAT data (fl3 ) and has been recently nicely confirmed by the PAMELA 
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observations ()82l ) of the positron fraction 

P+ <T) tot 

t = Z-S± fi 47) 

Different astrophysical contributions to the positron fraction in the 10 GeV 
region have aheady been been explored in (|l4l). More accurate and energy 
extended data will shed light on the effective presence of a bump in the 
positron absolute flux and on its physical interpretation. A possible expla- 



nation relies on DM species annihilating in the galactic halo (|2ll: l64l). This 



hypothesis has been strongly revived (|86l : l83l) with the PAMELA observation 



of what seems a clear positron excess above 10 GeV. 

Such an interpretation, though very exciting, is at some point limited by 
the uncertainties in the halo structure and in the cosmic ray propagation 
modeling. Equations ()1.28p and ()1.42p can be combined to yield the positron 
flux generated by the WIMP annihilations taking place within the Milky 
Way diffusive halo 

$°+^^(0, e = EjEo) = ^ X 4 X / des g{es) /dm (Ad) • (1-48) 
The information related to particle physics has been factored out in 

The energy distribution g{es) describes the positron spectrum at the source 
and depends on the details of the WIMP annihilation mechanism. The halo 
integral /dm is the convolution of the reduced positron propagator G with 
the square of the DM galactic density 

/dm(Ad)= / d3x5G'(xo^xs;AD) l^^l' . (1.50) 

In figure 11.91 the positron fraction ()1.47p is presented as a function of the 
positron energy E. The total positron flux at the Earth 

^tot ^ ^DM ^ ^sec (^ 5-^^ 

encompasses the annihilation signal and a background component for which 
the results of (figl ) as parameterized by (21) have been used - see the brown 
thin solid lines. The mass of the DM species is 100 GeV and a NEW pro- 
file has been assumed. The observations featured in the various panels are 
indications of a possible excess of the positron fraction for energies above 
10 GeV. Those measurements may be compared to the red thick solid curves 
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Fig. 1.9. The positron fraction e'^/{e + e+) is plotted versus the positron energy 
E for a 100 GeV WIMP in the case of a NFW BM profile. The four panels refer to 
different annihilation final states : direct e+e^ production (top left), bb (top right), 
W^W~ (bottom left) and t+t~ (bottom right). In each panel, the brown thin 
solid line stands for the positron background borrowed from (ITol) and parameter- 
ized by (21!). The red thick solid curve refers to the total positron flux where the 
signal is propagated with the best-fit choice of the astrophysical parameters, i.e. the 
configuration JVIED of table 11.21 The yellow area features the total uncertainty band 
arising from cosmic ray propagation. The different models found in ((SSl ) to be com- 
patible with the B/C ratio all yield a positron fraction which is enclosed inside this 
yellow strip. Experimental data from HEAT jll), AMS (0; [tI^AFRICE ^ 
and MASS (jH) are also presented for comparison. Figure from (jSSh. 



that correspond to the MED configuration of table 11.21 In order to get a 
reasonable agreement between the DM predictions and the data, the anni- 
hilation signal has been boosted by an energy-independent factor ranging 
from 10 to 50 as indicated in each panel. As is clear in the upper left panel, 
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the case of direct production offers a very good agreement with the positron 
excess. Notice how weh ah the data points he within the yehow uncertainty 
band. A boost factor of 10 is enough to obtain an excehent agreement 
between the measurements and the median flux. A smaher value would 
be required for a flux at the upper envelope of the uncertainty strip. The 
W~^W~ and t~^t~ channels may also reproduce reasonably well the obser- 
vations, especially once the uncertainty arising from cosmic ray propagation 
throughout the Milky Way diffusive halo is taken into account. They need 
though larger boost factors of the order of 30 to 40. On the contrary, softer 
production channels, like the bh case, are unable to match the features of 
the positron bump. 



1.6 Conclusions 

Crucial information on the astronomical dark matter pervading the Galaxy 
will be provided by crossing the positron f|82i ) and antiproton (jsil ) obser- 
vations. An excess seems to be present in the positron flux but does not 
show up in the antiproton signal. This may lead to strong constraints on 
the properties of the DM species as mentioned for instance in BH. Pure 
leptonic models are a priori favored since even a small branching ratio to 
quarks or gauge bosons would translate into an overproduction of antipro- 
tons which is not seen. Some attention needs to be paid though insofar as 
antiprotons could be rehabilitated depending on how cosmic ray transport 
is modeled. The forthcoming measurements will be crucial to ascertain the 
presence of WIMPs in outer space and to study their nature. 
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